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Abstract 

In this paper we elaborate on the idea [Lohse et al., Phys. Rev. Lett. 
78, 1359-1362 (1997)] that (single) sonoluminescing air bubbles rectify 
argon. The reason for the rectification is that nitrogen and oxygen 
dissociate and their reaction products dissolve in water. We give further 
experimental and theoretical evidence and extend the theory to other 
gas mixtures. We show that in the absence of chemical reactions (e.g., 
for inert gas mixtures) gas accumulation in strongly acoustically driven 
bubbles can also occur. 
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I. INTRODUCTION 



Sonoluminescence (SL) has long been known to be very sensitive to the used 
gas This effect is even more pronounced for single bubble sonoluminescence 

(SBSL), a phenomenon in which a single gas bubble is driven by a strong acoustic 
field and can emit short light pulses for hours [^0. Detailed experiments by the 
Putterman group at UCLA ||-|16| revealed that inert gas as part of the operating gas 
is essential. In figure ^ we show the experimental SL intensity from a SL bubble in 
water as a function of the percentage of noble gas mixed with nitrogen, taken from 



ref. |TT|. Pure nitrogen bubbles show hardly any SL emission; the optimum is around 



= 0.01 = 1%, i.e., the amount of argon contained in air. 




FIG. 1. SL intensity (normalized to air) from a SL bubble in water as a function of the 
percentage (mole fraction) of noble gas mixed with nitrogen. The gas mixture was dissolved 
in water at a pressure head of ISOmmHg, i.e., Poo/Po ^ 0.2. The data are taken from figure 



1 of Hiller et al. 



Besides the noble gas percentage C,i the other experimentally controllable param- 
eters in SBSL experiments are the forcing pressure amplitude Pa of the forcing sound 
field 

P{t) = Pa cos Ut, (1) 

the total gas pressure overhead Poo, and the ambient pressure Pq which is normally 
latm. The frequency 00/211 which is in a range between 20kHz and 40kHz is fixed 
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as it has to be adopted to the size of the resonator ("Crum cell", |T3)- Note that 
the ambient radius Rq of the bubble (i.e., the bubble radius at ambient standard 
conditions Pq = latm and at the temperature of the water) is not a free parameter, 
but the system adjusts Rq itself by diffusional processes. 

There are two types of SBSL: stable and unstable SBSL O, 



13. In unstable SBSL 



the phase of light emission (measured relative to the phase of the forcing sound field) 
drifts on a diffusional timescale of seconds, interrupted by sudden breakdowns. The 
same is true for the maximal bubble size per cycle and for the light intensity. An 
example for the bubble dynamics and light intensity for unstable SBSL is shown in 
figure 6 of ref . 
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In ref. [0 we quantitatively accounted for unstable SBSL as 
light emission from a bubble growing by rectified diffusion |]5|, p^,p!8|-p2| and thereby 
running into a shape instability [^, 0, where it pinches off a microbubble and 

starts over. The pinch off leads to a recoil of the bubble which makes it appear 



"jiggling" or 
visualized 



'dancing" 



Meanwhile, the pinched off microbubble could be 



p6| . The second type of SBSL is stable SBSL, distinguished by a constant 



phase and intensity of the light pulses, repeating for hours with remarkable precision 

mm- 



In ref. we calculated phase diagrams of SL bubbles in the ambient radius vs 
forcing pressure and gas pressure vs forcing pressure phase spaces. These diagrams 
are based on the Rayleigh-Plesset equation for the bubble radius R{t), a similar type 
of approximation for shape distortions, and the advection diffusion equation. We 
call this approach the Rayleigh-Plesset SL bubble approach. Shape stabilities and 
diffusional stabilities are considered. The phase diagrams quantitatively agree with 
Barber et al.'s [|1^ and Lofstedt et al.'s [|1^] experiments for argon (and other inert 
gas) bubbles, but not for air bubbles. 

For argon bubbles stable SBSL is only possible in a small window of tiny gas 
concentration, see figure ||. For Pa = L3atm this window is between / Pq = 0.002 
and 0.004 in very good agreement with Lofstedt et al.'s experimental data. Repeating 
the calculation of such phase diagrams for air bubbles essentially gives the same result, 
however, experimentally stable SBSL in air bubbles is found at about one hundred 
times larger gas pressure overhead [^, 13,2^. This discrepancy between air and argon 
bubbles was first pointed out by Lofstedt et al. |]T3[ who hypothesized an "as yet 
unidentified mass ejection mechanism" in air bubbles which "is the key to SL in a 
single bubble". 



II. AIR BUBBLES VS ARGON BUBBLES 



In refs. |]28|,|29[| we have suggested that this mechanism is chemical. The occurrence 
of chemical reactions has in fact led to the discovery of multibubble sonoluminescence 
(MBSL) ^0|-|32[| : Frenzel and Schultes were stimulated to look for lumines- 
cence as the formation of hydrogen peroxide in aqueous fluids subjected to sound 
had been observed before. Later, Schultes and Gohr Il33|l found that also nitric and 
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FIG. 2. Phase diagram for pure argon bubbles in the / Pq versus Pa/Po parameter 
space. Stable SL is only possible in a very small window of argon concentration. The 
experimental data points included for comparison refer to observed stable SL (filled symbols) 
or stable non-SL bubbles (open symbols) from refs. (diamonds) and [^2| (circles) and 
show good agreement with the theory. Note that only those data can be included for which 
Pa, Poo and S^i are experimentally known. 



nitrous acids were produced. The reason is that the high temperatures generated by 
the bubble collapse are beyond the dissociation temperature of oxygen and nitrogen 
(~ 9000K @]), leading to the formation of O and N radicals which react with the H 
and O radicals formed from the dissociation of water vapor. Rearrangement of the 
radicals will lead to the formation of NO, OH, and NH, which eventually dissolve in 
water to form H2O2, HNO2, and HNO3, among other products. 



Based on fits of SBSL spectra I10|,H1 and hydrodynamic calculations [T^,|3^, it is 



assumed that internal bubble temperatures in SBSL are even higher than in MBSL. 
Therefore, the same reactions as in MBSL will occur. The reaction products (NO2, 
NO, . . .) are pressed into the surrounding liquid, and are not recollected during the 
next bubble cycle, since their solubility in water is enormous. These chemical pro- 
cesses deprive the gas in the bubble of its reactive components. Small amounts of 
N2 and O2 that diffuse into the bubble during the expansion react and their dissoci- 
ation products are expelled back into the surrounding liquids at the bubble collapse. 
The only gases that can remain within a SBSL bubble over many bubble cycles are 
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those which at high temperatures do not react with the hquid vapor, i.e., inert gases. 
Hence, when air is dissolved in water, a strongly forced bubble is almost completely 
filled with argon. This argon rectification happens in SBSL but not in MBSL because 
it requires long time bubble stability. 

This argument immediately suggests that the partial pressure of argon 

= ilPoo (2) 

determines bubble stability, not the total pressure p^o- Indeed, if we include Holt 



and Gaitan's experimental data for air bubbles |22| on inert gas - nitrogen mixtures 



and only consider the inert gas partial pressure as the relevant quantity for diffusive 



stability, excellent agreement with the theoretical phase diagrams |]T8[ is found. 

The phase diagram shows that for Pa = 1.3atm argon bubbles exist between 
0.002 < p^ / Pq < 0.004. For pure argon dissolved in water, p-^ = p^. For air bubbles, 
however, the partial pressure of argon is only p"^ = O.Olpoo which requires 0.2 < 
Poo/Pq < 0.4 for stable SL with "air" bubbles, in good agreement with experiment. 
Since this amount of degassing is easily achieved, air with its 1% argon is a particularly 
friendly gas for SL experiments. The theory suggests that for an argon ratio of 
C,i ~ 0.0033 at Pa = 1.3atm there is stable SL between 0.6 < Poo/Po < 1-2, so that 
degassing is not required. Also, the window of stability is even wider than for air. The 
major problem in experimentally achieving SBSL without degassing is spontaneous 
cavitation, provoked by impurities in the liquid. These must be eliminated for the 
experiment to work. Another way to obtain stable SL without degassing would be to 
slightly increase the ambient pressure Pq so that the ratio p'^ /Pq is in the required 
window. Indeed, with such a kind of experiment stable SBSL without degassing could 
recently be achieved [p^ . 



Before we proceed with a quantitative analysis, we will give further support for the 
nitrogen dissociation hypothesis [p8|,p9[ from comparison with various experimental 
results on gas mixtures. 

• Transition towards SL: The transition towards SL with increasing forcing pres- 
sure Pa is shown in figure ^ for both argon and air bubbles. For pure argon 
bubbles the transition to SL is very smooth. For air, however, one can observe 
a breakdown of the bubble radius at about 1.1 atm, signaling that the dissocia- 
tion threshold of N2 is achieved. Before the transition the bubble is filled with 
a mixture of nitrogen, oxygen, and argon, and the ambient radius is determined 
by the combination of all three gases. Beyond the dissociation threshold, only 
argon is left in the bubble. - The transition from the no SL regime in figure 
to the SL regime can also be seen in figure 0, where the SL light intensity 
is plotted as a function of for fixed Pa (we assume Pa = 1.3atm) and fixed 
Poo/Pq ~ 0.20. According to our theory we expect SL for p^ / Po > 0.002 or 

> ^^75^ = 0.01 (3) 

Poo/Pq 
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FIG. 3. Transition in the maximal radius towards the SL regime for argon bubbles 
(dots) and for air bubbles (squares). Bubbles around Pa = l.latm start to glow. Only for 
air bubbles a breakdown in the radius is seen near that onset of SL, signaling the threshold 
for nitrogen dissociation. The data for argon are taken from figure 4 of Hiller et al.'s work 
|TT|| , the total gas concentration is about poo = 150mmHg. This means /Pq = 0.2 and 
according to figure the bubble should be in the unstable SL regime in agreement with the 



observations of ref. [11|. The data for air are taken from figure 2 of Barber et al.'s work |36|. 
The gas saturation is about 10%, corresponding to /Pq = 0.01 • 0.1 = 0.001. According 
to figure ^ we have stable SL around Pa ~ 1.4 atm which again is in agreement with the 
experiment reported in ref. [pq]. 



in pretty good agreement with figure || where we indeed observe that strong 
SL is "switched on" at about that concentration. Our theory also predicts that 
near the switch on we always have stable SL, whereas for larger unstable SL 
develops. The window of stable SL is only in between 0.002 < / Pq < 0.004, 
so that beyond = 0.004/poo we expect unstable SL, a prediction which should 
be verified. 

• Hysteresis: Above we have seen that by decreasing the percentage of inert 
gas down to ~ 0.0033 at Pa = 1.3atm, we obtain a very wide window of 
stable SL around Poo/Po = 1- When the percentage of inert gas is even 
lower, the argon partial pressure / Pq = ^iVoo/Pq can be pushed below the 
stable SL regime in figure ^ even for Poo/Po = 1- An example of this has been 
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FIG. 4. The transition (in Rmax) to SL for a bubble filled with an initial 0.1% xenon in 
nitrogen gas mixture at a partial pressure of ISOmmHg. Only the bubble at Pa = 1.4atm 
emits light. According to the dissociation hypothesis this corresponds to / Pq = 0.0002 
for a pure sonoluminescing xenon bubble. From the phase diagram figure ^ which is (with 
tiny corrections) also valid for xenon we conclude that bubbles driven at Pa = 1.3 atm 
dissolve for these low concentrations, whereas bubbles at Pa = 1.4 atm show stable SL, just 
as seen here. The data of this figure are taken from figure 11 of Lofstedt et al. |15|. 



observed in figure 11 of Lofstedt et al. [|I^, see figure ^ of the present work, 
which studies SBSL of xenon doped nitrogen bubbles (which behave like argon 
doped nitrogen bubbles) in water at S,i = 0.001 and Poo/Po = 0.2, corresponding 
to p^/Pq = 0.0002, provided that the gas temperature is high enough to exceed 
the dissociation temperature of nitrogen. This is the case at Pa = 1.3atm and 
according to figure ^no stable SL is possible. However, for larger Pa = 1.4atm 
stable SL becomes possible again (as the window of stable SL moves down), 
just as observed in experiment. Moreover, the system shows hysteresis: the SL 
state at Pa = 1.4atm can be reached by continuously increasing the forcing and 
thus continuously substituting N2 with xenon. On the other hand, decreasing 
the forcing pressure below the stability threshold for pure xenon bubbles leads 
to the dissolution of the bubble. 
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Unstable SL: An example for unstable SL is shown in figure 6 of ref. For 
that figure we have = 0.05 and Poo/Po = 0.20. Thus, from equation (H) we 
have /Pq = 0.01 and according to the phase diagram figure ^ we are well in 
the unstable SL regime, just in agreement with the observations. 

Large Pa bubble dissolution islands: Further support comes from the recent work 



of Holt and Gaitan p2[ who measured detailed phase diagrams as a function 
of the ambient radius Rq and the forcing pressure Pa for different air pressure 
overhead Poo/Po- Their central result is that at Poo/Po = 0.2 there is a relatively 
large forcing pressure Pa ~ 1.2— 1.3atm regime where bubbles dissolve, see figure 



Such dissolution islands do not exist within theories of rectified diffusion 15 



T8|-p0|, which predict the bubbles to grow in that regime. However, considering 
the nitrogen dissociation at large Pa, one realizes that the observed dissolution 
islands are a direct consequence: chemical reactions deplete the bubbles from 
air and only argon is left. As the relevant partial pressure p^ / Pq = 0.002 is so 
low, bubbles must indeed shrink in that regime. 

Isotope scrambling: Also experiments with hydrogen gas support the dissocia- 
tion hypothesis. Hiller et al. [0 analyze SL in H2 and D2 gas bubbles, both 



in normal and in heavy water. Since the gas dynamics inside the bubble deter- 
mines the strength of the light emission, the SL intensity curves should group 
according to the gas content. However the four (H2 in H2O, H2 in D2O, D2 in 
H2O, and D2 in D2O) experiments group according to the surrounding liquids, 
cf. figure 2 of ||12[|. This suggests the following scenario: Both the gas and the 
liquid vapor in the bubble dissociate to some extent during the (hot) compres- 
sion phase and recombine later on during expansion. Even if there is only a 
minor amount of this isotope "scrambling" , after thousands of cycles the gas in 
solution around the bubble would contain the same hydrogen isotope as the bulk 
liquid. This kind of "scrambling" is well known in MBSL PSI. The acoustic 



resonator theory of SBSL ||3^ stipulates that the light intensity decreases with 
increasing acoustic transmission, which scales with the ratio of the gas density 
to the liquid density. Since D2 is heavier than H2, SBSL in heavy water should 
be dimmer than in normal water, as observed by Hiller and Putterman [|l^ . 



III. MODELING THERMAL EFFECTS WITHIN RAYLEIGH-PLESSET 

BUBBLE DYNAMICS 

We now proceed to a quantitative calculation of phase diagrams for gas mixtures. 
The dynamics of the bubble radius R{t) is well described by the Rayleigh-Plesset 
equation §,|D[. 

RR + ^if = - {p{R, t) - P{t) - Po) 
^ Pi 
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FIG. 5. This figure is adopted from Holt and Gaitan's measurement figure Ic of ref. 
It is the phase diagram in the Ro-Pa parameter space for air at Poo/Po = 0.20; the 
driving frequency is 20.6kIIz. Arrows indicate whether the bubbles grow or shrink. Three 
equilibrium curves A, B, and C can be recognized. In between curves B and C there is a 
"dissolution island" . The shaded area shows the shape stable parameter domain. Note that 
in the SL region the onset of shape instabilities is at Rq ~ 7fim rather than at Rq ~ 5/im 
as theoretically calculated in ref. |]l^ . One contribution to this deviation is that the driving 
frequency chosen by Holt and Gaitan is smaller than the w/27r = 26.5kHz used by Barber 
et al. and throughout the calculations presented here; another one is the oversimplified 
model of thermal effects we use. 



R d , ^ R 2a 
+ —-Mt)-Au--— (4) 
pici at R piR 



with a van der Waals pressure 
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FIG. 6. Forcing pressure (in atm) with amplitude Pa = 1.3atm, bubble radius R{t) 
(in fim) with the typical collapse ("Rayleigh collapse") after the maximum, Mach number 
M{t) = R{t)/cg{t) (where Cg is the speed of sound in the gas), and temperature 0(t) (in 
K) as function of time for a bubble with Rq = 2^m (top to bottom). 



Typical parameters for an argon bubble in water at room temperature [|T3| 



are 



the surface tension a = 0.073kg /s^, the water viscosity u = 10~^m^/s, density 
pi = 1000kg /m^, the polytropic exponent 7 = 1, and speed of sound in water 
Q = 14:81m/ s. For all calculations in this work we picked the same driving frequency 
as in the SL experiments performed on argon bubbles |jl3[, namely ci;/27r = 26.5kHz, 



corresponding to a period T = 27t/uj = 38/is. Finally, h = i?o/8-86 is the hard core 
van der Waals radius for argon bubbles . Typical time series for the bubble radius 



R{t) resulting from equation (^) for given forcing P{t) are shown in fig. 

Thermal conduction effects can approximately be taken into account by simply 
putting the polytropic exponent 7 = 1 in equation (^. The reason is that the bubbles 
employed in SL experiments are so tiny and the oscillation period T ^ 38 fis so long 
that the gas in the bubble equilibrates with the water temperature G;. A quantitative 
analysis was carried out by Prosperetti (see also the review [Q) who calculated 
how the polytropic exponent 7 depends on the (thermal) Peclet number Pe = RIu/k, 
see figure 1 of ref. For experimental support of this analysis, see Crum [|^. The 
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Peclet number gives the ratio between the bubble length scale Rq and the thermal 
diffusion length k/uj. The thermal diffusivity n for argon is k ^ 2 • lQ~^rn? / s, which 
yields Pe 0.2 for Rq = 5fim and according to figure 1 of ref. |42], the poly tropic 
exponent 7 = 1. 

The RP equation obviously contains much smaller time scales than u^^. One 
could therefore argue that these smaller timescales may enter into the calculation of 
Pe, so that the frequency u should be replaced by \R\/R. This estimate leads to an 
instantaneous Peclet number 

which can become as large as as 10^ at the Rayleigh collapse. According to figure 
1 of ref. showing 7(Pe) or to Prosperetti, Crum, and Commander's fit of that 
curve [0 this implies 7 ~ 5/3 for argon at the time of the collapse. However, since 
Pe{t) ^ 1 only holds in very small time intervals ~ Ins, the global dynamics are 
hardly affected by setting the effective polytropic exponent 7 = 1 uniformly in time. 
Note that with 7 = 1 equation should not be thought of as an equation of state but 
rather as a process equation parametrizing the isothermal conditions at the bubble 
wall, induced by the large heat capacity of water. The choice of 7 = 1 is confirmed by 
the full numerical simulations of Vuong and Szeri |H5| and by the analysis of Kamath 



et al. Note that, as a consequence, there are heat fluxes back and forth across 
the bubble wall. 

Another confirmation for putting 7 = 1 for the bubble dynamics comes from a 



recent work by Yasui who modeled the heat flow through the bubble wall and 



also included water evaporation and condensation. Consideration of these effects [^| 



fits the experimentally measured R{t) curves [ll3|, p!4|, p6| pretty well However, 
the essence of the dynamics is the same as with simply putting 7 = 1 and the phase 
diagrams we are going to calculate in this paper are hardly affected by neglecting the 
details of the thermal effects. 

Even if thermal effects are simply treated by putting 7 = 1 for most of the time 
during the bubble oscillation, they are essential at the final stage of the collapse and 
thus for estimates of the temperature achieved in the bubble. The present understand- 
ing H, 44] is that the only way to reliably predict the temperatures in the bubble is a 



full numerical simulation of the gas dynamical equations inside the bubble. Clearly, 
for strong enough forcing shock waves will bounce back and forth inside the bubble 
and the spatial temperature distribution in the bubble will be highly inhomogeneous 



Such full numerical calculation is not within the spirit of the Rayleigh-Plesset 
SL bubble approach We therefore come back to Prosperetti's [Hj] model and 



estimate also in the case of nonlinear bubble oscillations the temperature within that 
model. I.e., we calculate the instantaneous Peclet number (Rf), determine the resulting 



effective polytropic exponent jeff from figure 1 of |^ and assume a van der Waals 
behavior of the internal gas according to 
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e(t) = e, 



7<ij?-l 



(7) 



where G/ is the temperature of the hquid. This, as any model which neglects spatial 
inhomogeneities in the gas, is a very crude model for the temperature and will only be 
able to give orders of magnitude. The physical statements resulting from it, however, 
do not depend on details of this model and so we feel that such an approach is justified. 
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FIG. 7. The minimal radius, the maximal radius, the maximal effective polytropic 
exponent 7e//, and the maximal temperature Qmax are shown for two different ambient 
radii Rq = 2fim (solid) and Rq = 5/im (dashed) as a function of the forcing pressure 
amplitude Pa- For smaller bubbles the crossover from the harmonically oscillating regime 
for small Pa to the bouncing bubble regime for large Pa is more abrupt. 

It can be seen from figures 0c and |^d that this temperature model gives the correct 
trends. In figure 0c we show the maximal •jeff per cycle as a function of Pa- We also 
give the minimal and the maximal bubble radius in figs. ^ and b. For small Pa 
the maximal effective polytropic exponent 7e// is close to one and jumps towards 
5/3 at the transition from the sinusoidal to the bouncing bubble R{t) dynamics. 
This transition has been analyzed in detail in ref. |^0|- The temperature 0(t) from 
equation (0) is shown in figure |^d. It strongly peaks at the coUapse; for all other 
times it essentiahy equals the liquid temperature 9/ (= 293K in this example). The 
heating at the collapse can be enormous. 
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FIG. 8. Maximal temperature in the bubble as a function of Rq and Pa within our 
model calculation. 

What are the temperatures achieved in this approach? In figures ^ and ^ we 
plot the maximal temperature within our model as a function of Rq and Pa- In the 
parameter regime of interest for SL, values as high as 20000 — 60000-ft" are achieved. 
These values agree order of magnitude wise with the independent results of Bernstein 



et al. |3^,|51| who extracted the temperature from the spectral shape and obtained 



30000 to 60000K, depending on the gas contents. 

IV. PHASE DIAGRAMS FOR AIR BUBBLES 

Now consider a bubble in water containing a mixture of a reactive gas (taken to 
be N2) and an inert gas Ar. The total number of moles of gas in the bubble is 

^^ot = ^^ = N^,+N^^, (8) 

where 0o = 273K is the normal temperature and G = 8. 31 A3 J /{molK) is the gas 
constant. The argon ratio in the bubble is 

& = ^. (9) 



that of nitrogen 



i-& = ^. (10) 



If c^^{r,t) and c'^^{r,t) are the concentration fields of Ar and N2 in the liquid, re- 
spectively, the rate of change of the number of moles of N2 and Ar in the bubble is 
given by 
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I\Ar = (-LJ-J 



r 



Nn2 = ^^iV2 exp -— . (12) 

Here, -Dat-, Djq^, jiAr and /iTVa are the respective diffusion constants and molecular 



masses. The concentration fields obey a mass advection diffusion equation |T8|,pO 
whose boundary conditions are set by the external concentrations 

-TO 

(Henry's law) and by the partial gas pressures p"(t) in the bubble 



c"{R{t),t) = c\ 



■ 



Pn 



a = Ar, N2. The solubilities for nitrogen and argon are different, Cq = Q.Q61kg/m 
and Cq^ = 0.020kg /m^ |^2|. The diffusion constants are approximately the same 



D^r = Dfyf^ = 2- 10~^m^/s. The second term in (|T2p represents the bubble's nitrogen 
loss by chemical reaction. The reaction rate will depend on the temperature G(t) in 
the bubble. For simplicity, we assume that the reactions follow an Arrhenius law, 
with empirical parameters appropriate for nitrogen dissociation (ref. |0): A = 6 ■ 
10^^(6o/6)^-^(po//^Af2)(-Ro/-R)^cm^/(mo/-s) giving the timescale of the reaction; 6* = 
IISOOOK ^1| is the activation temperature and po the equilibrium gas density. This 
reaction law is rather crude, as it neglects backward reactions as well as the kinetics 
of the expulsion of reaction products; however, it is sufficient for this demonstrative 
calculation. 

We can straightforwardly extend the adiabatic approximation of the slow diffu- 
sional dynamics [|^,[l^,|2^ (i.e., separation of times scales) to equations (pT]), ([T^). 



The only requirement is that the involved chemical reactions are fast compared to 
diffusional processes which definitely is the case. The result of the adiabatic approx- 
imation is that the change per cycle is given by 



ANm2 ^tiDn^Cq 



'^{Pi'-^b{p),) (13) 

7^j|c^ _ _ J _ ^^^^p (-e7e))o (14) 

with the weighted time averages 

/o R''{t)dt 

For further simplification in the numerical calculations to come we in addition employ 
the saddle point approximation of the integral I defined by 



14 



as / ~ 1/Rmax, as already done in [jT5|. 

The last term in eq. (H) will only contribute at the collapse when G(t) is large. 
The chemical reaction rate will depend on the detailed space and time dependence 
of the temperature in the bubble which we model by the polytropic law (0) - as 
explained in the previous section. 
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FIG. 9. Phase diagram for air at Poo/Po = 0.20 in the Rq - Pa space. The arrows denote 
whether the ambient radius grows or shrinks at this parameter value. Curve A denotes the 
equihbrium for an air bubble, on curve C the bubble only contains argon. The intermedi- 
ate curve B necessarily exists because of the topology of the diagram, and represents an 
additional stable equilibrium. The thin line shows when the nitrogen dissociation threshold 
9000K is reached. 



With these approximations equilibrium points AA^^^ = AA^at^ = in the two 
dimensional space {N^r, N^^), or equivalently in the space {^b,Ro), can easily be 
calculated from the RP dynamics via equations (]T3|), (|T3|). As in ref. |T^, only 
averages of type (|15D have to be determined for each parameter pair (i?o, Pa)', there- 
fore, phase diagrams can be calculated solely from the Rayleigh-Plesset equation (^). 
The equilibrium radii Rq in the Rq — Pa plane for air = 0.01) at Poo/Po = 0.20 
are shown in figure H. For small forcing the temperatures are not high enough to 
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initiate chemical reactions, so that the equihbrium curve corresponds to the predic- 



tion of ref. [|l^] for this gas concentration. This equihbrium is unstable: The bubble 
either shrinks or grows by rectified diffusion. As pointed out above the growing bub- 
ble eventually runs into a shape instability where microbubbles pinch off and make 
the bubble dance because of the recoil [|I^]. In the opposite limit of high forcing 
(curve C), the reactions burn off all the N2, so that the bubble contains pure argon; 
this equilibrium corresponds to the (stable) equilibrium at the argon partial pressure 
/Po = O.Olpoo/A = 0.002. 

Figure ^ displays a regime of shrinking bubbles at high forcing pressures (left of 
curve C) and an adjacent region of growing bubbles (right of curve A). This neces- 
sitates the existence of an additional equilibrium at intermediate forcing pressures, 
curve B in figure ^, for which growth by rectified diffusion and mass loss by reactions 
balance. This additional equilibrium occurs close to the point of nitrogen dissocia- 
tion, and turns out to be stable; the argon fraction for this equilibrium is slightly 
larger than the fraction in the liquid (for not too strong forcing). 
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FIG. 10. The top curve shows the equilibrium argon ratio Q m a Rq = 5/xm bubble as 
a function of the forcing pressure Pa- The bottom curve presents the corresponding argon 
and nitrogen loss/gain per cycle, normalized to the total amount of moles. 



In figure [T^ we present the net loss/gain per cycle for both argon and nitrogen 
for a Rq = 5fim bubble. At small Pa < l.lSatm the second term in equation ( |14|) 
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is not important as the gas does not become hot enough. The bubble diffusively 
shrinks (for Pa < 1.045atm) or grows because of rectified diffusion; the argon ratio 
is constant at = = 0.01. In between the shrinking and the growing regime is 
the unstable equilibrium point A. Around Pa = 1.17atm the heating is sufficient to 
lead to some nitrogen dissociation and at Pa = l.lSatm the nitrogen loss through 
dissociation balances its gain through rectified diffusion. At that forcing pressure 
also argon is in equilibrium as the argon ratio has increased up to ~ 0.1 so that 
Poo — C {p)4- have reached the stable equilibrium point B. At larger forcing 
pressure Pa the bubble is in the dissolution island. For increasing Pa the nitrogen 
loss becomes less because less and less nitrogen is left in the bubble. The argon loss, 
on the other hand, diminishes because we approach the stable diffusive equilibrium 
point C at P(j ^ 1.4atm. Beyond this forcing pressure, the bubble grows again by 
rectified argon diffusion; essentially no nitrogen is left in the bubble (1 — pa 6 • 10~^ 
at Pa = 1.40atm). Note that because of the separation of time scales the values on 
refer to an average value over the whole cycle. When the bubble is close to its 
maximum, the nitrogen concentration will be slightly larger, right at the collapse the 
argon concentration will be slightly larger. 




FIG. 11. The fraction of argon in the bubble as a function of Pa and Rq. 



How robust is this picture, i.e., how does it depend on details of temperature 
dynamics and the chemical reaction mechanism which are both poorly understood? 
The large Rq parts of branches A and C do not depend at all on these details. The 



reason becomes clear from figure [TT| in which we show the composition equilibrium 

(17) 



as a function of Rq and Pa- It is given by 
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Plugging (|T3|) and (JT^ into (0) gives the quadratic equation for the equihbrium 
argon ratio 



a + e)a' + {/3-e)e,+S = (18) 



with the coefficients 



a = ^Mi(^-^:^), (19) 

3Geo ey/o' ^^^^ 

^ = - (22) 





J/i 



Ar 



From the two solutions only the one with < .^^ < 1 has physical meaning. The 
interpretation of figure |ll| is straightforward: Weakly forced bubbles have ~ 6, 
thus Poo/Pq = 0.20 is relevant for stability. Strongly forced bubbles have ~ 1 ^ 
thus p^' / Pq = 0.002 is the relevant quantity. The transition between these regimes is 
abrupt, and occurs when the bubble temperature surpasses the dissociation temper- 
ature (~ 9000/^ for N2). Where exactly this happens depends on the model of the 
temperature and thus so does the equilibrium curve B in figure ^ 

However, the existence of this additional stable equilibrium is independent of the 
model details. It simply follows from topological reasons: To the right of the curve A 
the bubbles are growing, to the left of curve C they are shrinking, so in between there 
must be an equilibrium. Indeed, as mentioned already above and shown in figure 



Holt and Gaitan's ||22| recent detailed measurements of Ro — Pa phase diagrams found 
such a classically unexpected equilibrium (for the same air pressure Poo/Po = 0.20 as 
chosen here) and the connected large Pa "dissolution island" between curves B and 
C. 



Holt and Gaitan repeated the experimental phase space measurements for 



larger air pressures Poo/Po = 0.40 and Poo/Po = 0.50. The theoretical phase diagram 
for Poq/Pq = 0.50 and C,i = 0.01 (air) is shown in figure |12[ The island of dissolution 



between curves B and C is now smaller, as the partial argon pressure is p^ / Pq = 0.005 
and curve C is thus shifted to the left, whereas the nitrogen dissociation curve is hardly 
affected and consequently neither are curves A and B. In experiment we see the same 
qualitative behavior, but even stronger pronounced, see figure la of ref. [^. We 
expect the poor modeling of both the temperature dependence and the chemistry as 
the origin of the quantitative disagreement of the theoretical and experimental phase 
diagrams. 

For comparison we also calculated the phase diagram for Poo/Po = 0.20 and = 
0.025 and plotted it into the same figure O. For this situation the argon percentage 
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FIG. 12. The solid lines (curves A, B, C) show the same phase diagram for air as in 
figure ^, but now for Poo/Po = 0.50, i.e., /Pq = 0.005. For comparison, the open symbols 
(curves A', B', C) show the phase diagram for a nitrogen-argon mixture with an argon ratio 
of ^/ = 0.025 at Poo/Pq = 0.20, which beyond the nitrogen dissociation threshold results in 
the same partial argon pressure p^ /Pq = 0.005. Indeed, for large Pa the phase diagrams 
agree. Only curves A and A' below the nitrogen dissociation threshold differ. 

/ Pq = iiVoolPf) = 0.005 is the same as for Poo/Po = 0.50, C,i = 0.01. Indeed, the 
phase diagram demonstrates very convincingly that for large forcing pressure it is 
only the partial argon pressure which is the relevant parameter for diffusive stability: 
The equilibrium curves B and C for the two cases agree very well. 

V. PHASE DIAGRAMS FOR FURTHER NITROGEN- ARGON MIXTURES 

Our theoretical method can be applied to many other experimental situations. 
As a further example we calculate the corresponding phase diagram figure ^ to the 
experimental situation in figure 11 of [T^, see figure ^ here, where Poo/Po = 0.20 
and C,i = 0.001. The total gas concentration Poo/Po has the same value as in figure 
^. Therefore, the curves A and B are also essentially as in that figure, because they 
are determined by the nitrogen concentration. Curve C which is beyond the nitrogen 
dissociation threshold, however, is further to the right, reflecting the lower argon 
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FIG. 13. Rq - Pa phase diagram for nitrogen with a small amount = 0.001 of argon 
dissolved. The relative pressure overhead is Poo/Po = 0.20. These parameters correspond 
to the experimental situation in figure 11 of Lofstedt et al. jl5|, see figure § of the present 
work. 



concentration p^' / Pq = 0.0002. 

What happens if Pa is slowly increased? For small Pa = latm the bubble will 
be repelled from the unstable equilibrium A. It either dissolves or grows by rectified 
diffusion and runs into the shape instability where it pinches off a microbubble and 
the whole process starts over, as described in detail in ref. |]TB[. For Pa ~ l.lSatm the 
bubble does not run into the shape instability any more but in the attractive diffusive 
equilibrium curve B where it stays. On further increase of Pa it follows that curve and 
shrinks drastically. Exactly this sequence of events has experimentally been observed, 
see figure 11 of Lofstedt et al. [0. 

For larger Pa on first sight there seems to be a contradiction to the experimental 
observation figure 11 of |15|. Theoretically, the bubble should always stay on the stable 
equilibrium curve B. Experimentally, however, between Pa = l.Satm and 1.4atm no 
bubble seems to exist and at Pa = 1.4atm the ambient radius is much larger, similar 
to that of curve C rather than that of curve B. 

This variance may be resolved when we consider the unavoidable jitter in the 
forcing pressure amplitude Pa- At Pa = 1.35atm, the equilibrium curves A and B are 
so close that a small jitter APa (to smaller Pa) of less than 0.05atm (which is less than 
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the precision to which Pa can be determined; also the translational movements of the 
bubble and its distortion on the pressure field could contribute to a jitter.) suffices 
to make the bubble jump to the unstable domain below curve A where it dissolves. 
Thus, though the equilibrium B in that regime is stable from a mathematical point of 
view, from a physical point of view it is not and the bubble can dissolve. Support for 
this claim could come from an experiment in which the jitter is artificially increased: 
the gap in figure ^ where no bubble exists should widen. 

With the help of the phase diagram figure |13] we now better understand the above 
mentioned hysteresis. The stable part of branch C is reached by a fast increase of 
Pa- The bubble is then boosted into the growing regime below the stable branch 
of curve C. Subsequently, this stable branch C is reached through rectified diffusion. 
On the other hand, if Pa is decreased from Pa = 1.4atm, the bubble will first follow 
the stable equilibrium curve C and shrink. At the bifurcation point it falls off the 
equilibrium curve and shrinks towards the stable curve B. Because of the jitter in 
Pa curve B may not be able to trap the bubble as explained above and the bubble 
dissolves. Even a fast decrease of Pa towards much smaller values ~ l.latm does not 
save the bubble from dissolution as it only contains argon and the nitrogen needed for 
a stable equilibrium has been burned off in the first part of the cycle during which Pa 
was increased and recollecting nitrogen only happens on the slow diffusive timescale. 

For even lower argon concentration as e.g. in figure where we chose C,i = 0.0001 
and Poo/ Pq = 0.20 the situation is similar. The only difference is that curve C moves 
even further to the right and may not be reached in experiment because of the (short 
timescale) shape instabilities discussed in refs. [p!^, p5|] . 

Let us discuss the situation in the limiting case = 0, i.e., pure nitrogen. As 
explained above for = 0.001, on increasing Pa the bubble should first "dance" 
because of growth by rectified diffusion and pinching off of microbubbles and then 
it should be caught by the stable (but not necessarily sonoluminescing) equilibrium 
curve B. At this equilibrium bubble growth by rectified diffusion and nitrogen loss by 
dissociation again just balance. For large Pa ~ 1.3atm the bubble should dissolve by 
the same mechanism discussed above: Curves A and B are so close that a jitter in Pa 
leads to a jump from the stable equilibrium curve B to the dissolution domain below 
the unstable curve A. We suggest to experimentally look for the stable equilibrium 
curve B which is suggested by the nitrogen dissociation hypothesis. 

What experiments on pure nitrogen bubbles have been done up to now? To our 
knowledge only one experiment has been reported (Hiller et al. |ll[). In that experi- 
ment the nitrogen was 99.7% pure. Hiller et al. [|ll[| managed to measure extremely 
weakly sonoluminescing unstable bubbles which were hard to keep alive. This sup- 
ports the presented theory that in the strong forcing regime the high temperature 
destroys the nitrogen and thus the bubbles. However, the light intensity figure 5 of 
ref. [|ll| shows an oscillating pattern (on a timescale of seconds) which we do not 
understand. 
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FIG. 14. Phase diagram in Ro-Pa space for Poo/Po = 0.20 and a tiny argon concen- 
tration of = 0.0001. Only for very high forcing pressure diffusively stable bubbles are 
possible (but are likely to be unstable towards the surface shape instability). The bottom 
figure shows an enlargement of the lower right part of the top figure. 
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VI. INERT GAS MIXTURES 



It is interesting to note that the chemical instabihty of one component of the 
gas mixture is not a necessary requirement for the accumulation of one gas species 
in the bubble. A similar accumulation can be achieved for inert gas mixtures if 
the diffusion constants and the equilibrium concentrations of the two gases in the 
mixtures are different. As an example we consider a 1:1 mixture of helium and 
argon at Poo/Po = 0.20, i.e., = 0.50 for the argon ratio. The diffusion constants are 
Dav = 2- 10~^m^/s and Dhs = 5.8- 10~^m^/s |^| and the equilibrium concentrations 



are Cq^ = 61g/m^ and c^*^ = 0.78g/m^ The molecular masses are fi^r = 4:0 g/mol 
and iiHe = 4g/mol. The resulting equilibrium concentration in the bubble is shown 
in figure |l^. For small forcing helium accumulates in the bubble and for large forcing 
argon accumulates. The borderline between these two regimes is the unstable diffusive 



equilibrium curve calculated in |T8| and denoted "curve A" within the plots shown 
here. As no reactions occur, this is the only equilibrium: Below that curve bubbles 
diffusively shrink, above it they grow by rectified diffusion. From this behavior the 
helium or argon accumulation can immediately be qualitatively understood: The 
diffusive change of mass is proportional to the product of the material constants 
Dco/fi, see eq. (|T3|). Now DatCq^ / f^Ar = 3.1 ■ 10~^ mols^^m~^ and DhsCq^ / Hue = 
l.l-10~^mo/s~^m~^. Thus in the shrinking regime argon can escape faster and helium 
accumulates. Vice versa, in the growing regime the rectified diffusion of argon into the 
bubble wins and it accumulates. In the process of helium accumulation the helium 
equilibrium concentration 1 — ^^ may never be achieved as the bubble dissolves earlier. 
This problem does not occur for the process of argon accumulation, as the diffusively 
growing bubble runs into the shape instability and can survive after microbubble 
pinchoff [g. 

The maximal argon concentration achieved is only ^^,max ~ 0.73, less than the 
result Qmax = 1 of the previous section where the chemical instability of nitrogen 
caused its complete extinction. Also the transition in fig. |15| is less abrupt than in 
above case fig. |Tl| for which chemical reactions occurred. 

The maximal argon concentration can easily be calculated analytically from equa- 
tions ( pTsD - (1^) with the index N2 replaced by the index He. Obviously, e = as 
no reactions take place. From figure 7 of we realize that for large forcing (p)^ is 
small compared to p^ = p^ = 0.1. With this large forcing pressure approximation 
we obtain 



0.73 (23) 



in agreement with the numerical result figure [15 . 



23 



Pjatm 

FIG. 15. Argon equilibrium ratio in a bubble as a function of ambient radius Rq 
and forcing pressure Pa- The external gas is a 1:1 mixture of argon and helium at an 
overhead pressure ol poo/Po = 0.20. Only the difference in the diffusion constants and in 
the equilibrium concentrations between argon and helium leads to the argon accumulation 
in the bubble for large forcing. 

VII. CONCLUSION, PREDICTIONS, AND SUGGESTED EXPERIMENTS 



We conclude this paper with a summary of the nitrogen dissociation theory p8|,p9 



which we elaborated here in detail. It accounts for the dependence of SBSL on the 
percentage of inert gas in the bubble. It is based on a combination of principles from 
sonochemistry and hydrodynamic stability. The main result is that strongly forced air 
bubbles eventually consist of pure argon because at the high temperatures achieved 
in the bubble the nitrogen and oxygen molecules dissociate and react to water soluble 
gases. Consequently, it is the partial pressure of argon (or of any other inert gas) 
which is relevant for stability p9[. The dynamics of the bubble is described by the 
Rayleigh-Plesset equation; chemical reactions are assumed to obey an Arrhenius type 
law; back reactions are neglected. The bubble pressure is approximated by a van der 
Waals law. The slow diffusive processes can be treated in an adiabatic approximation. 
The theory contains no adjustable fit parameter. 

Finally, we make detailed suggestions of experiments to further test the nitrogen 
dissociation hypothesis. 

• Phase diagrams in the Poo/Po "^^s Pa/Po parameter space: The most basic pre- 
diction is that in the large forcing regime the partial pressure p^ /Po is the 
relevant parameter for obtaining stable SL (or the partial pressure of any other 
inert gas) . We suggest to map out the phase diagram figure |^ for argon-nitrogen 
mixtures at various pressure overheads Poo/Po and various argon ratios ^f. (i) 
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Only the partial argon pressure / Pq = ^iVoo/Pq should matter and (ii) the 
measured phase diagram should look like the theoretical one figure ||. 

Stable SL without degassing: In particular, stable SL will be possible without 
degassing, i.e., at Poo/Po = 1, if the argon ratio is properly adjusted. For Pa = 
1.3atm it must be ~ 0.0033. To test this prediction, extremely pure water has 
to be used to avoid spontaneous cavitation. An alternate way to achieve stable 
SL without degassing is to increase the ambient pressure correspondingly. 

Oxygen- argon mixtures: Oxygen has a dissociation temperature which is only 
slightly lower than that of nitrogen. Therefore, oxygen-argon mixtures should 
behave very similar to nitrogen-argon mixtures. In particular, at large forcing 
only argon is left in the bubble, so even the light intensities and spectra should 
be very similar. This prediction has recently been confirmed in experiment 



Other inert gas doped molecular gases: Apart from CO, nitrogen is the most 
stable molecular gas. Other molecular gases dissociate earlier. Therefore, mix- 
tures of them with inert gases will show a slight shift of the stable equilibrium 
curve B towards smaller forcing pressure. As the temperature increase with 
forcing pressure is very sudden (see figure |^), this shift will, however, be very 
small and may be undetectable. 

Stable nitrogen bubbles: Another very basic prediction of the nitrogen dissocia- 
tion theory is the existence of the equilibrium curve B in the phase diagrams of 
type 1^. For air Holt and Gaitan have observed it |2^. We predict that it exists 



for any reactive gas. In particular, we expect that pure nitrogen bubbles can 
be stable. 

Island of bubble dissolution for large Pa: Between the equilibrium curves B and 
C, bubbles dissolve. The size of this large Pa island of dissolution should only 
depend on the partial argon pressure p^ / Pq = C,iPoo/Po, hardly on the argon 
ratio C,i or the pressure overhead Poo/Po individually. For larger p^ / Pq the size 
of the dissolution regime should shrink, for smaller p^ / Pq it should grow. 

Radon doped SL bubbles: By using radon-nitrogen gas mixtures it may be pos- 
sible to detect the inert gas accumulation in the bubble because of the radioac- 
tivity of radon. 

Inert gas accumulation: Even for a mixture of nonreactive g He and Ar we 

in general have C,b ^i- In the large forcing pressure regime argon accumulates 
for this example. 

pH measurements for air bubbles: A very appealing test is to measure the 
concentration of the reaction products of the dissociated gases as a function of 
time, as already done in MBSL PS[. For nitrogen-argon mixtures nitrous acid 
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production would lead to a decrease in pH as already predicted in ref . . For 
an estimate of the production rate we assume that all the nitrogen that diffuses 
into the bubble during the bubble expansion is burned off at the collapse. This 



amount is estimated in ref. |15] as 



N2 



(24) 



per cycle. With typical values of Rmax = 10i?o for the maximal radius, Rq = 
5/im, Dn^ = 2 - IQ-^m^/s, ^ 0.20c^^ Cq' = 0.02kg /m^, and T = 38/is 
one obtains AN^^ ^ 3 • 10^^^ mol per cycle or ~ 3 ■ 10~^° mol of N2 per hour 
converted to reaction products. The consequence is a small but detectable pH 
decrease: In a 100m/ flask with pure water (pH=7) there are initially 10~* mol 
H"*" ions. Assume that all ejected N atoms eventually form either HNO2 or 
HNO3, then 6 ■ 10~^° mol H'*' ions are produced. This means a concentration 
increase to (10~^ + 6 ■ 10~^)mo/// or a pH value of 6.975 after an hour. 

• pH measurements for bubbles of other gas mixtures: As a blind test one should 
also measure the pH as a function of time for pure argon bubbles. No pH change 
should occur. Neither do we expect a pH change for oxygen-argon mixtures as 
the main reaction product of the dissociated oxygen should be H2O2. 
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